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1. INTRODUCTION 

With the increase of the population around the world, electrical energy management becomes very 
important, especially in the effort to enhance the stability of the power system [1]. From the perspective of 
angle stability, one of the best ways to dampen the swing and improve stability is the excitation control using 
Flexible Alternating Current Transmission Systems (FACTS) devices. Static Var Compensator (SVC) is one 
of the widely used FACTS devices for power system applications. In addition to being used to control 
reactive power and supporting bus voltage, SVC is also capable of enhancing angle stability in the system 
[2, 3]. Over the past few decades, there are various techniques that have been proposed for controlling SVC 
to improve the damping of synchronous engines. Among the methods proposed are proportional-integral (PI) 
controllers [4, 5], proportional-integral-derivative (PID) controllers [6-9] and lead lag (LL) controllers 
[10, 11]. In this paper, a combination of SVC with PID controller (SVC-PID) has been chosen to improve 
oscillation stability in power system. 

With the improvement of AI technology, optimization techniques have become an option in 
managing control parameters and solving other issues in power system stability. Among them are 
Evolutionary Programming (EP) [12-14] and Artificial Immune System (AIS) [15-17]. Both EP and AIS 
share almost identical features, whereas EP mimicking biological evolutionary processes, AIS works on the 
principle of human immune system. Although both algorithms are concise and do not use complex process 
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streams, EP and AIS are not capable to provide optimal solutions and take longer computation time. In this 
paper, PSO has been proposed to implement the optimization process. PSO work with a simple optimization 
concept, resulting an easy implementation, manage to provide best optimal solutions and consumed shorter 
time to converge. 

In searching for the optimal results, the selection of indicator used as the objective function in the 
optimization process is one of the important criteria. The indicator should be sensitive enough to evaluate the 
system, so that the optimal value of parameters can be optimized. One index is chosen as an objective 
function for conventional optimization techniques. Although it is successfully enhanced the damping 
capability, single index only improves certain property in eigenvalues. From eigenvalues point of view, 
usage of minimum damping ratio, Gin as objective function can shift eigenvalue region towards the real axis 
of phase plan by the decreasing value of imaginary part of the eigenvalue @ [18, 19]. Meanwhile, maximum 
damping factor, Oia, shows better result compared to ¢,,i, in capability of pushing the set of eigenvalues 
further left of the imaginary (jm) axis of the phase plan [20, 21]. To encompass both advantages, one multi- 
objective function (MO) is proposed. 

This paper presents the impact of MO in PSO based angle stability improvement technique for 
SMIB system. The study involved an effective approach from objective function point of view, to determine 
the parameters of SVC damping controller in solving angle stability problems. PID controller has been 
chosen for the damping controller attached to the SVC system and its fixed-gains are determined using PSO 
optimization technique. The goal is to stabilize the system in minimum time. Along with designating MO, 
previously used objective functions; Gin aNd Onax are selected for comparison purpose. 


2. SINGLE MACHINE CONNECTED TO AN INFINITE BUS (SMIB) SYSTEM 

Single machine connected to an infinite bus (SMIB) system is an equivalent system with one 
synchronous machine connected to one big bus with an infinite load. To enhance the damping capability, 
SMIB is equipped with a flexible AC transmission (FACTS) device called static var compensator (SVC). 
Widely used nowadays, SVC has very fast control response time, capable to reduce the need for maintenance 
of the absence of rotating parts and manage to perform individual phase control [2]. In this paper, 
PID controller is installed together with SVC to attempt better damping performance. 

Figure | shows the Phillips-Heffron block diagram model for the SMIB system that equipped with 
SVC and PID controller (SVC-PID) [6]. SVC-PID can be segregated to SVC and PID element. SVC element 
involve of circuit constant, Ky and time constant, Ty. The PID element is obtained from proportional gain, Kp, 
integral gain, K; and derivative gain, Kp of PID controller. Before carrying out Ao to SMIB system, the value 
of SVC-PID have to be limited in specified upper and lower limits. 
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Figure |. Phillips-Heffron block diagram model of the SMIB system with SVC-PID controller 


Here, H is the inertia constant, K, is the damping torque coefficient, Tp and Kr are the time constant 
and circuit constant of the exciter oscillation system, respectively. The value of wo is equal to 27fo. K), Kz, K3, 
T;, K4 Ks and Ks are constants associated with the operating real and reactive load, impedance ratio, 
electrical torque, excitation level, rotor speed and rotor angle in the generator. Detail calculation of the 
system can be found in [4]. 
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3. MATHEMATICAL MODELLING FOR SMIB SYSTEM WITH SVC-PID 
Based on the Phillips-Heffron model of the SMIB system equipped with SVC-PID as shown in 
Figure 1, the following equations are applied: 
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As shown in (1)-(6) can be encompassed into a matrix form as followed: 
Xsvc = Asve * Xsvc + Bsvc *U (7) 
where 
Xsyvc = [4@, 46 AE, Ave Ao Ap] (8) 


Comprehensive calculation of the system in Figure | can be found in [4]. Parameter values for the 
generator, transmission line, exciter and SVC controllers of the SMIB system are selected based on 
parameters in [6]. 


4. MULTI-OBJECTIVE FUNCTIONS FOR ANGLE STABILITY IMPROVEMENT 

To elevate the damping capability of power systems, one multi-objective function (MO) derived 
based on combination of minimum damping ratio, ¢,;, and maximum damping factor, o,,,, 1s proposed. 
The Gnin 1S an indicator that gives the minimum value among the damping ratio of every eigenvalues in the 
SMIB system, based on the real and imaginary part of every eigenvalues at the loading condition. With the 
optimization of ¢,nin, the system poles are consistently pushed further left of the imaginary (j@) axis. Also, the 
decreasing value of imaginary part of the eigenvalue @ will shift eigenvalue region towards the real axis. 
Effects of these combinations can be shown as a wedge-shaped sector on the phase plan [18]. The regions of 
eigenvalues on the phase plan, before and after optimization process are showed in Figure 2(a). 
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Figure 2. Comparison of eigenvalue areas on the phase plan 
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The Oma, calculates the maximum value of damping factor o of every eigenvalue in the 
electromechanical modes of oscillation. The optimization of Gq, will result the system poles to be shifted 
further left of the jw axis [20]. The regions of eigenvalues on the phase plan; before and after optimization 
process are showed in Figure 2(b). By implementing the proposed MO, area of eigenvalues will shift further 
left of the imaginary axis and at the same time towards real axis. The area of eigenvalues on the phase plan 
for comparison case between system with MO approach and without MO approach is shown in Figure 2(c). 
The area bounded by MO effect can be shown as a trapezoid-shaped sector on the phase plan. To combine 
both nin aNd Onax Into one fitness equation, an objective function is designated as follows: 


Juo = % Simin + 2° Omax = %* min (-«/ of + oF) + a@2*Max(Gmax),0 S a4,a2 <1 (9) 


a, and a2 are constants; attached to Gnin and Oma, fespectively to tune the percentage of both indicators. o; and 
«@; are respectively the real and imaginary part of the i” eigenvalue. The combination of the objective 
function can be formulated as: 


Maximize (Jno) (10) 
Jno Tefers to multi-objective function (MO). Following limitations must be complied for SVC-PID system: 
Ke < Kp < Kee Ky" < K, < Kr Ky < Kp < Ki al 1) 


To estimate the maximum value of J), the following process was implemented: 
a) Calculate eigenvalues using (7) 
b) Evaluate J,,. using (9). 
c) Vary the value of selected parameters and repeat step (i) and (ii) until the difference between the 
maximum and minimum value of objective function, Jing, and Jinin, respectively reaches 0.001 or number 
of iterations reach 100. 


5. PARTICLE SWARM OPTIMIZATION (PSO) TECHNIQUE 

Particle Swarm Optimization (PSO) technique is developed based on the movement of flock of birds 
flocking and school of fish in searching for food. PSO tries to adopt the principle of individual intelligence as 
well as the influence of collective behavior of the herd. This algorithm involved initialization, velocity and 
positions update, fitness calculation and best position update. The flow chart which represents PSO algorithm 
is illustrated in Figure 3. Completed descriptions of PSO algorithm process can be found in [16]. 
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Figure 3. Flowchart for the particle swarm optimization (PSO) algorithm 















































Here, Kp, K; and Kp are optimized by PSO method. The parameters involved in the proposed 
approach is tabulated in Table 1. 


Table 1. The Parameters for PSO Algorithm 


List of Parameters C1, = C2 = 0.5, @max = 0.09, Omin = 0.04 
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6. RESULTS AND ANALYSIS 

In this paper, two cases were implemented to search the most suitable ratio of Gi, and Oma, that 
compress in one multi-objective function (MO). Both cases were conducted using single machine connected 
to an infinite bus (SMIB) system in MATLAB environment. The simulation used SVC-PID system with 
Particle Swarm Optimization (PSO) as the computational intelligence technique. All values of Kp, K; and Kp 
parameters are randomly selected. The simulations are simulated in two different loading conditions as 
tabulated in Table 2. For loading conditions which are less or more than the values simulated in Table 2, 
they are considered to produce the similar results. 


Table 2. Three Different Loading Conditions for Case | and Case 2 








Case 1 2 
Loading Active power, P (p.u.) 0.35 0.75 
condition Reactive power, Q (p.u.) -0.2 0.5 





Emin ANd Omax have different value characteristics. The value of €,i, is in the range of [0,1]. On the 
contrary, there are large possibility of value for Oia, (from -0o to +00). To acquire the benefits of both 
indicators, the merger of the two indexes is essential. For this to happen, the o,,,, is reduced to maximum 10% 
of actual value, so one index will not monopolize other index in the new indicators. 

In searching for the most suitable value of Gin and Omax, twelve different MO combinations have 
been introduced based on the (9). Comparison of constants a; and a2 for all twelve different MO 
combinations for all cases are tabulated in Table 3. 


Table 3. Comparisons of Ratio Constants a, and a, for all Twelve MO Combinations 








MO Combinations Baneconsians MO Combinations Reto constants 
Qi A2 Qy Ae 

Zt 1 0 MO-6 1 0.06 
MO-1 1 0.01 MO-7 1 0.07 
MO-2 1 0.02 MO-8 | 0.08 
MO-3 1 0.03 MO-9 1 0.09 
MO-4 1 0.04 MO-10 1 0.10 
MO-5 1 0.05 Sg 0 0.10 


To analyze the results for all three cases, all responses are divided into two groups: Group | and 
Group 2. Group | consist of result for Zt, MO-1, MO-2, MO-3, MO-4, MO-5 and MO-6, while Group 2 
cover result for MO-6, MO-7, MO-8, MO-9, MO-10 and Sg. Figure 4(a) and 4(b) show the response of speed 
deviation of Group | and Group 2 for Case 1, respectively. 
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(a) Group 1 (MO-1~MO-5 and Zt) (b) Group 2 (MO-6~MO- 10 and Sg) 


Figure 4. Speed response in various combination for Case 1 


From Figure 4(a) and 4(b), it shows that all combinations are completely damped within 1.6 seconds 
of total computation time, which indicate that all twelve approaches demonstrate stable condition. From the 
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speed deviation response, the oscillation patterns are decreased from MO-1 to MO-6 approaches, gradually, 
and increased again from MO-7 to MO-10. From the comparison of responses, it is discovered that MO-6 
gives the fastest damping responses and less oscillation, followed by MO-7. Zt approach shows almost the 
same result with MO-1. On the other hand, Sg gives the biggest oscillation rate, which makes it the worst 
among all the twelve combinations. 

The regions of eigenvalues location in complex s-plane for Group 1 and Group 2 respectively are 
shown in Figure 5(a) and 5(b). Regarding eigenvalues, the one located at the left-hand side of s-plane, which 
justifies that all the twelve system approaches are stable. Among those twelve combinations, the eigenvalues’ 
region based on MO-6 approach is located far to the left side of s-plane and closer to the real-axis. This result 
revealed that the MO-6 combination is the most prominent approach from the rest. Meanwhile, Zt and Sg 
approaches are among the worst combinations. In spite of positioned at the left-hand side area, the region of 
eigenvalues for Zt is at the nearest position to the right-hand side of s-plane compared to the other eleven 
combinations. For Sg, the eigenvalues distributions are located further away from the real-axis. 


$ 
On @ 0% 





Cc a : 
joie: * wx Sos « 
Sms Cg * ace * 

nr = 5 Q Se Cf 

2 e 

: ® 

ES g 
ca ey 

“10 g “10 * 

a 7 ” Real Patt (p.u.) . i > " 7 Real Part (p..) 
(a) Group 1 (MO-1~MO-5 and Zt) (b) Group 2 (MO-6~MO-10 and Sg) 


Figure 5. Location of eigenvalues for Case | 


The comparisons of SVC-PID parameters based on twelve multi-objective (MO) combinations for 
Case | tabulated in Table 4. The optimized values of Kp for all twelve combinations are in narrow range, 
between 0.5039 and 0.5693. The same pattern of results is also obtained for K; (-0.0414~-0.0005) and Kp 
(0.1461~0.5611). This shows that all twelve combinations tuned the SMIB systems almost in the same 
manner, with slightly different value. This also indicates that the tuning process of damping parameters is 
very dependent on the selection of the objective function. In this critical process, only the best combination 
will calculate the most suitable SVC-PID parameters and gives the most improvement in stability condition 
of SMIB system. 


Table 4. Comparisons of SVC-PID Parameters Based on MO Combinations for Case | 








st SVC-PID Parameters meus SVC-PID Parameters 
MO Combinations Kp K Kp MO Combinations Kp K, Kp 

Zt 0.5106 -0.0212 0.1461 MO-6 0.5282 -0.0384 0.3324 
MO-1 0.5362 -0.0174 0.2129 MO-7 0.5192 -0.0332 0.3517 
MO-2 0.5693 -0.0183 0.2721 MO-8 0.5262 -0.0210 0.4131 
MO-3 0.5255 -0.0064 0.3528 MO-9 0.5510 -0.0070 0.4622 
MO-4 0.5242 -0.0046 0.4031 MO-10 0.5039 -0.0005 0.5611 
MO-5 0.5649 -0.0111 0.3877 Sg 0.5315 -0.0414 0.2626 





The speed deviation response of Group | and Group 2 for Case 2 are shown in Figure 6(a) and 
Figure 6(b), respectively. As those in the previous case, all twelve MO combinations are capable to elevate 
the stability condition of SMIB system. In this case, MO-3, MO-4, MO-5, MO-6 and MO-7 are the top five 
achievers in giving the low oscillation rate in speed deviation response. In other hand, Zt and Sg approaches 
are among the worst, resulting bigger oscillation compare to other combinations. This shows that the two-in- 
one objective function is very convincing in giving better damping capability, rather than a single objective 
function. 
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Figure 6. Speed response in various combination for Case 2 


The regions of eigenvalues location in complex s-plane for Group 1 and Group 2 for Case 2 are 
shown in Figure 7(a) and 7(b). In this case, the regions of eigenvalues for all twelve combinations are 
positioned at the left-hand side of s-plane. Among those twelve, MO-6 is succeeded to shift its eigenvalues 
region to the farthest on the left side of s-plane and most close to real-axis. This justifies that the combination 
of MO-6 is the most suitable approach for Case 2. On the other hand, the approach of Sg gives good result by 
positioning the eigenvalues coordinate to the far left of s-plane. But, as two of the eigenvalues are situated at 
the farthest from real-axis, it is making Sg not so efficient to improve the stability condition of the system. 
Meanwhile, one of eigenvalues of Zt method is positioned at the nearest to origin but at the right-hand side of 
s-plane, resulting Zt as the worst objective function approach. 
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Figure 7. Location of eigenvalues for Case 2 


Table 5 tabulates the comparisons of SVC-PID parameters based on MO combinations for Case 2. 


Table 5. Comparisons of SVC-PID Parameters Based on MO Combinations for Case 2 








MO Combinations _ sala lai Zz MO Combinations a e 
Zt 0.6281 -0.3887 0.0815 MO-6 0.5232 -0.3028 0.1285 
MO-1 0.6124 -0.2900 0.1276 MO-7 0.5773 -0.2021 0.1775 
MO-2 0.5732 -0.1988 0.1458 MO-8 0.5870 -0.1744 0.1983 
MO-3 0.4527 -0.4492 0.0932 MO-9 0.6021 -0.1191 0.2705 
MO-4 0.5130 -0.4235 0.0988 MO-10 0.6025 -0.1535 0.2185 
MO-5 0.5822 -0.3457 0.1193 Sg 0.6235 -0.1077 0.2794 








Oscillation stability enhancement using multi-objective swarm based technique for... (N. A. M. Kamari) 


638 o ISSN: 2502-4752 


Two cases presented have proved that the combination of maximum damping ratio €,i, and 
maximum damping factor Oy, , aS multi-objective function is succeeded to enhance the optimization 
capability to become more accurate in tuning SVC-PID parameters. Almost all the combinations which have 
been presented demonstrate better damping results as compared to conventional single objective functions 
such as Gnin ANd O;ax. AMong those ten combinations, the MO-6 (a; = 1, a= 0.06) approach shows the most 
prominent results based on two responses; speed deviation and eigenvalues region. Using this approach, 
PSO optimization technique constantly identify the optimal values of Kp, K; and Kp to decrease the 
oscillation rate of speed deviation. In complex s-plane result, the MO-6 approach is also capable to shift the 
eigenvalues position to the most left of left-hand side of imaginary-axis and to the nearest to real-axis. 


7. CONCLUSION 

From this paper, it can be concluded that introduction of multi-objective (MO) manage to elevate 
the stability quality of power system. The application of minimum damping ratio €,,;, as objective function of 
an optimization process will shift the eigenvalues’ region further to left side of complex plane, while the 
usage of maximum damping factor oj, will bring eigenvalues position towards the real axis of s-plane. 
Merger of Gin ANd Oma, Index in the certain ratio as objective function will bring both advantages into the 
selected optimization process. From Cases | and 2, both results indicate that combination of Gin and Omax 
approach at a, = | and a2 = 0.06 elevate the performance of PSO computation technique, given the best 
damping capability of SVC-PID controller among all approaches. 
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